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JOINT FORMATION OF QSOS AND SPHEROIDS 

G.L. Granato^, L. Silva^, G. De Zotti^, and L. Danese'^ 
RESUMEN 

El resumen sera traducido al espanol por los editores. In view of the extensive evidence of tight 
inter-relationships between spheroidal galaxies (and galactic bulges) with massive black holes hosted at their 
centers, a consistent model must deal jointly with the evolution of the two components. We describe one 
such model, which successfully accounts for the local luminosity function of spheroidal galaxies, for their 
photometric and chemical properties, for deep galaxy counts in different wavebands, including those in the 
(sub)-mm region which proved to be critical for current semi-analytic models stemming from the standard 
hierarchical clustering picture, for clustering properties of SCUBA galaxies, of EROs, and of LBGs, as well as 
for the local mass function of massive black holes and for quasar evolution. Predictions that can be tested by 
surveys carried out by SIRTF are presented. 

ABSTRACT 

In view of the extensive evidence of tight inter-relationships between spheroidal galaxies (and galactic bulges) 
with massive black holes hosted at their centers, a consistent model must deal jointly with the evolution 
of the two components. We describe one such model, which successfully accounts for the local luminosity 
function of spheroidal galaxies, for their photometric and chemical properties, for deep galaxy counts in different 
wavebands, including those in the (sub)-mm region which proved to be critical for current semi-analytic models 
stemming from the standard hierarchical clustering picture, for clustering properties of SCUBA galaxies, of 
EROs, and of LBGs, as well as for the local mass function of massive black holes and for quasar evolution. 
Predictions that can be tested by surveys carried out by SIRTF are presented. 
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1. INTRODUCTION 

The hierarchical clustering model with a scale in- 
variant spectrum of density perturbations in a Cold 
Dark Matter (CDM) dominated universe has proven 
to be remarkably successful in matching the observed 
large-scale structure as well as a broad variety of 
properties of galaxies of the different morphological 
types (Granato et al. 2000 and references therein). 
However, serious shortcomings of this scenario have 
also become evident in recent years. 

The critical point can be traced back to the rel- 
atively large amount of power on small scales pre- 
dicted by this model which would imply far more 
dwarf galaxies or substructure clumps within galac- 
tic and cluster mass halos than are observed (the 
so-called "small-scale crisis"), unless star formation 
in small objects is strongly suppressed (or the small 
scale power is reduced by modifying the standard 
model). At the other extreme of the galaxy mass 
function we have another strong discrepancy with 
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model predictions, that we might call "the massive 
galaxy crisis" . Even the best semi-analytic mod- 
els hinging upon the standard picture for struc- 
ture formation in the framework of the hierarchical 
clustering paradigm, are stubbornly unable to ac- 
count for the (sub)-mm (SCUBA and MAMBO, Fig- 
ure 1) counts of galaxies, most of which are probably 
massive objects undergoing a very intense starburst 
(with star formation rates ~ 1000 Mq yr"-'^) at z > 2. 
Recent optical data confirm that most massive ellip- 
ticals were already in place and (almost) passively 
evolving up to z ~ 1-1.5. These data are more con- 
sistent with the traditional "monolithic" approach 
whereby giant ellipticals formed most of their stars 
in a single gigantic starburst at substantial red- 
shifts, and underwent essentially passive evolution 
thereafter. In the canonical hierarchical clustering 
paradigm the smallest objects collapse first and most 
star formation occurs, at relatively low rates, within 
relatively small proto-galaxies, that later merged to 
form larger galaxies. Thus the expected number of 
galaxies with very intense star formation is far less 
than detected in SCUBA and MAMBO surveys and 
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Fig. 1. Integral source counts at 850 /xm and 1.4 mm predicted by the model by Granato et al. (2001) compared with 
observations. The dotted, dashed, dot-dashed lines show the contributions of starburst, spiral, and forming elliptical 
galaxies respectively. The long-dashed (at 1.4 mm only) gives the counts of radio sources (Toffolatti et al. 1998). 



the surface density of massive evolved ellipticals at 
^ > 1 is also smaller than observed. The "mono- 
lithic" approach, however, is inadequate to the ex- 
tent that it cannot be fitted in a consistent scenario 
for structure formation from primordial density fluc- 
tuations. 

2. RELATIONSHIPS BETWEEN QUASAR AND 
GALAXY EVOLUTION 

The above difficulties, affecting even the best cur- 
rent recipes, may indicate that new ingredients need 
to be taken into account. A key new ingredient may 
be the mutual feedback between formation and evo- 
lution of spheroidal galaxies and of active nuclei re- 
siding at their centers. In this framework, Granato 
et al. (2001) elaborated the following scheme. 

Feed-back effects, from siipernova explosions and 
from active nuclei (note that supernova feedback 
alone falls short of solving the dearth of dwarf galax- 
ies, hut photo-ionization by the UV background re- 
ionizing the inter-galactic medium (IGM) could do 
the job) delay the collapse of baryons in smaller 
clumps while large ellipticals form their stars as soon 
as their potential wells are in place; the canonical hi- 
erarchical CDM scheme - small clumps collapse first 
- is therefore reversed for baryons. Large spheroidal 
galaxies therefore undergo a phase of high (sub)-mm 
luminosity. 

At the same time, the central black-hole (BH) 
grows by accretion and the quasar luminosity in- 
creases; when it reaches a high enough value, its 
action (ionization and heating of the gas) stops the 



star formation and eventually expels the residual gas. 
This explains the observed correlation between BH 
and host spheroidal masses. The same mechanism 
distributes in the IGM a substantial fraction of met- 
als and may pre- heat the IGM. The onset of quasar 
activity (and the corresponding end of star forma- 
tion) occurs earlier for more massive objects. The 
duration of the starburst increases with decreasing 
mass from ~ 0.5 to ~ 2 Gyr. This implies that the 
star formation activity of the most massive galaxies 
quickly declines for z <3, i.e. that the redshift dis- 
tribution of SCUBA/MAMBO galaxies should peak 
at z > 3, as quasars reach their maximum luminos- 
ity (at z ~ 2.5). This explains why very luminous 
quasars are more easily detected at (sub)-mm wave- 
lengths for z > 2.5. 

A "quasar phase" follows, lasting 10^-10^ yrs, 
and a long phase of passive evolution of galaxies 
ensues, with their colors becoming rapidly very red 
[Extremely Red Object (ERO) phase]. Intermediate- 
and low-mass spheroids have lower Star Formation 
Rates (SFRs) and less extreme optical depths. They 
show up as Lyman-Brcak Galaxies (LBGs). 

Therefore, in this scenario, large ellipticals evolve 
essentially as in the "monolithic" scenario, yet in 
the framework of the standard hierarchical cluster- 
ing picture. Many aspects and implications of this 
compound scheme have been addressed by our group 
in a series of papers (Granato et al. 2001, Maglioc- 
chetti et al. 2001, Perrotta et al. 2002, Romano et 
al. 2002). Here we only summarize how the scenario 
compare with sub-mm counts. 
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Fig. 2. Integral counts and z distribution for a flux limit of 30|UJy at 24 /jm by the model by Granato et al. (2001). 
Sybols for counts as in Figure 1. The three-dots/dashed line shows the total counts of ellipticals, including also those 
where the star formation has ended. In the right-hand panel, the solid and the dotted lines show the z distributions of 
ellipticals during and after the star formation phase, respectively. 



2.1. Counts at (sub)-mm wavelengths 

The (sub)-inm counts are expected to be very 
steep because of the combined effect of the strong 
cosmological evolution of dust emission in spheroidal 
galaxies and of the strongly negative K-correction 
(the dust emission spectrum steeply rises with in- 
creasing frequency). The model by Granato ct al. 
(2001) has extreme properties in this respect: above 
several mJy its 850 /zm counts reflect the high-mass 
exponential decline of the mass function of dark ha- 
los. In this model, SCUBA/MAMBO galaxies corre- 
spond to the phase when massive spheroids formed 
most of their stars at z > 2.5; such objects essentially 
disappear at lower redshifts. On the contrary, the 
counts predicted by alternative models (which are 
essentially phenomenological) while steep, still have 
a power law shape, and the redshift distribution has 
an extensive low- 2; tail. As illustrated by Figure 1, 
the recent relatively large area surveys are indeed 
suggestive of an exponential decline of the 850 /um 
counts above several mJy. Further evidence in this 
direction comes from MAMBO surveys at 1.2 mm. 



2.2. Predietions for SIRTF surveys 

SIRTF surveys have the potential of providing 
further tests of the model. In particular the 24 /xm 
survey to be carried out as a part of the GOODS 
(http://www.stsci.edu/science/goods) Legacy Sci- 
ence project should reach a flux limit of 100 /zJy. Ac- 
cording to the model, about 50% of detected galaxies 
should be spheroidal galaxies forming their stars at 
z > 2. About 400-600 such objects are expected 
over an area of 0.1 square degree (Figure 2). Their 
redshift distribution is predicted to peak at z slightly 
above 2, with a significant tail extending up to 2; ~ 3. 
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